Spin-flip transitions induced by time-dependent electric flelds in surfaces with strong 

spin-orbit interaction 
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We present a comprehensive theoretical investigation of the light absorption rate at the 
Pb/Ge(lll)-/3^/3 x V3i?30° surface with strong spin-orbit coupling. Our calculations show that 
electron spin-flip transitions cause as much as 6% of the total light absorption, representing one order 
of magnitude enhancement over Rashba-like systems. Thus, it is demonstrated that a substantial 
part of the light irradiating this nominally non-magnetic surface is attenuated in spin flip processes. 
Remarkably, the spin-flip transition probability is structured in well defined hot spots within the 
Brillouin zone where the electron spin experiences a sudden 90° rotation. This mechanism offers 
the possibility of an experimental approach to the spin-orbit phenomena by optical means. 

PACS numbers: 71.18.-|-y, 73.20.r, 73.25.-^i 



Understanding electron spin transport and spin relax- 
ation in quasi-two-dimensional (2D) systems is of a capi- 
tal importance due to both, the fundamental reasons and 
the potential technological applications. The spin-orbit 
(SO) interaction is the most prominent relativistic effect 
leading to the fascinating phenomena recently observed 
in 2D systems, such as the quantum spin Hall effect ■ 
An experimentally accessible spin degree of freedom of- 
fers new route for the emergent field of spintronics, where 
the main features of charge dynamics are strongly influ- 
enced by the spin- related effects 0]. The technical 
possibility of spin manipulation and control by means of 
an applied bias voltage is strongly supported by recent 
investigations on a variety of semiconducting alloy sam- 
ples HQ. 

However, a strong SO coupling cannot be achieved in 
conventional semiconductors where the spin-splitting of 
the conduction electrons is limited to a few meV at most. 
As opposed, the relativistic effects completely dominate 
the electronic structure of many heavy-element surface 
materials and overlayers Q. The reason resides on the 
breaking of the inversion symmetry and the associated 
gradient of the effective one-electron potential introduced 
at the interface. These effects lead to extraordinarily 
large (100-500 meV) spin-splittings among the so called 
Shockley-type surface states, as it was first observed in 
the free electron-like Au(lll) noble metal surface @. 
This system is regarded as a prototype of the standard 
Rashba model [13, [HI ■ Considering the general form of 
the SO interaction in the non-relativistic limit. 
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the Rashba SO coupling is recovered by taking the gradi- 



ent of the effective one-electron potential V{y) as a con- 
stant and completely surface perpendicular. This model 
produces an entirely isotropic result, with a simple lin- 
ear spin-splitting of the two spin sub-bands with chiral 
spin-polarization. However, many heavy-element adlayer 
or even clean surface materials, such as Bi(ll O) Il2i . 
lxlH/W(110) 0L_A.u/Si(557) Q, Bi/Si(lll) 
or Tl/Si(lll) [Utill exhibit large anisotropic SO inter- 
action, inducing complex spin textures that considerably 
deviate from the free electron-like picture of the Rashba 
model. 

In this Letter, we investigate electrically induced spin- 
fiip excitations on the Pb/Ge(lll)-/3A/3 x A/3i?30° sur- 
face (\/3Pb/Ge(lll)), considering the full spinor struc- 
ture of the electron states within ab-initio density func- 
tional theory. The goal is to understand and quantify the 
striking mechanism leading to spin-fiip transitions out of 
a purely electric perturbation in a non-magnetic material. 
This system presents two well defined spin-spit surface 
states crossing the Fermi level, while the bulk substrate 
remains semiconducting. Thus, we face a problem involv- 
ing a completely spin-polarized 2D electron gas which is 
essentially decoupled from the bulk, i.e., an optimum sce- 
nario for studying surface spin-flip transitions [20| . Note- 
worthy, the strength of the SO interaction associated to 
the Pb atoms is two orders of magnitude larger (r^ 300 
mcV) than in semiconductor quantum wells. Moreover, 
the Pb overlayer includes strong anisotropic gradients in 
the surface ionic potential, giving rise to a fast variation 
of the noncollinear spin-polarization in momentum space 
which critically enhances the spin-flip transition proba- 
bility. 



The spin manipulation in 2D systems l21| is accessible 
via electric-dipole spin resonance mechanism, which cou- 
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FIG. 1: (color online) Top view of the V3Pb/Ge(lll) sur- 
face [2^. The small (gray) spheres symbolize the Ge sub- 
strate layers, whereas the big (green) spheres represent the 
Pb surface monolayer. The solid (red) parallelogram indi- 
cates surface unit cell. 



pies the spin-dependent electron velocity to externally 



applied fields [23. At surfaces, the difi'erent spin-spit 



sub-bands are connected by inter-band transitions that 
flip the electron spin [2^ . The detailed analysis of these 
processes provides valuable information about the spin- 
dynamics of the system, thus making the first- principles 
approach imperative. Simplified tight-binding Rashba- 
like models, although very useful for understanding the 
basic physics 2J], appear not to be realistic enough to 
account for material specific details which determine the 
electromagnetic response. 

Let us consider the interaction of an electron with 
an external time-dependent electric field of frequency w. 
Within the electric dipolc approximation valid for small 
momentum transfer q 0, the interaction Hamiltonian 
is given by 2^ [23 1 
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where Acxt(i) — A.QCOSiot is the external vector poten- 
tial associated to the electric field Ecxt {t) = Eq sin cot 
with Eg = Aqw/c, and v is the electron velocity opera- 
tor. In systems with SO coupling, besides the canonical 
contribution p/m, the velocity operator includes an ad- 
ditional spin-dependent term which emerges as the main 
responsible for the spin-fiip transitions [23|. 

Within first order perturbation theory, the transition 
rate associated to spin-fiip excitations between spin-spit 
surface states 5* and S' due to photon absorption (term 
Aoe~*'^ */2) can be derived from the Fermi's golden rule. 
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The integral is taken over the surface Brillouin zone 
(SBZ), fi^i^w) and e^^k represent the Fermi-Dirac distri- 
bution and a surface state eigenvalue, respectively, and 
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FIG. 2: (color online) Electron band structure of the 
\/3Pb/Ge(lll) surface. The scalar and fully relativistic bands 
are represented by dashed (red) and solid thick (blue) lines, 
respectively. The continuous background is the bulk band 
projection. The fully relativistic metallic surface states are 
labeled as S and S' . 



AIss' (k) is the interband matrix element 
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Above, \E'ik(i") denotes the single-particle Bloch spinor 
wave function associated to a surface state. 

Due to the inherent phase indeterminacy of the Bloch 
states in k-space, the matrix elements of the velocity op- 
erator require a special treatment pif . Following the 



approach presented in Ref. [27|, we expressed Mss'(k) 



in Eq. [4] in terms of the so called maximally localized 
Wannier functions [2^ . In this way, the matrix elements 
entering Eq. [3] are maximally smooth and thus suit- 
able for any interpolation procedure within the Brillouin 
zone. The Wannier functions were generated considering 
the entire structure of ab-initio spinor wave functions 
obtained within the noncoUinear DFT formalism (29I - 
3l| . The exchange-correlation energy was a ppr oximated 
by the standard LDA-PZ parametrization [32| and the 
2x2 norm-conserving fully relativistic pseudopotential 
approach [si, 33|. The ground state self-consistent cycle 
was performed considering the usual Monckhorst-Pack 
mesh corresponding to a 27 x 27 grid. We employed a 
very fine 200 x 200 mesh considering the standard Wan- 
nier interpolation procedure for all the ingredients enter- 
ing Eq. El [13, IS] in order to reliably account for the 
details close to the Fermi level jiTj . 

The V3Pb/Ge(lll) surface (Fig. [J) was simulated by 
a repeated slab technique containing 14 Ge layers. The 
Pb monolayer was included only in one side of the slab, 
while the other (bare) Ge(lll) surface was covered by a 
hydrogen adlayer in order to saturate the dangling bonds. 
We also analyzed the Au(lll) metal surface considering 
a 22 Au layer slab. In both systems, a full geometry op- 
timization was performed until all atomic forces exerted 
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FIG. 3: (color online) (a) and (b) Momentum dependent spin-polarization structure associated to the S and S' surface states, 
respectively. Arrows represent the in-plane spin-polarization component, whereas the background indicates the surface perpen- 
dicular component of the magnetization, mz(k). The Fermi surface of each state is indicated by solid (blue) lines, (c) Spin-flip 
transition probability associated to the S and S' surface states for R-circularly polarized light. 



on individual atoms were negligibly small (< 10 ^ Ry 
a.u.-i). 

Fig. [5] shows the calculated band structure of the 
Pb/Ge(lll)-\/3 surface. While the scalar relativistic 
calculation shows a single spin-degenerate surface band 
crossing the Fermi level, fully relativistic calculations 
present two spin-polarized surface bands labeled as S 
and S' . The SO interaction has a huge impact on the 
electron structure close to the Fermi level, such that this 
term cannot be treated pcrturbatively. Its contribution 
is even more important than some non-relativistic DFT 
terms such as the exchange-correlation energy. In this 
context, the SO interaction completely determines the 
metallic character of the S and S' surface states. These 
exist as surface states only outside the area close to the 
r point, where S and S' become resonances entering the 
bulk projection (continuum in Fig. [2]). Outside this re- 
gion, the splitting is in overall of the order of 100 meV, 
reaching a maximum of 300 meV near the high symme- 
try point M. The calculated Fermi wave vectors along 
the high symmetry direction FM, kp ~ 0.41 and 
kp ~ 0.37 A~^, are in very good agreement with re- 
cent ARPES experiments reporting kp = 0.40 and 
k§! = 0.36 A-i 

A characteristic feature emerging from the SO interac- 
tion at surfaces is the momentum-dependent spin polar- 
ization: 

m„(k) = J *;k(r)<T*nk(r)d3r, (5) 

where n is the band index. In Figs. [3^ and[3j3 we present 
the calculated spin-polarization for the S and S' surface 
states in the entire Brillouin zone. These figures demon- 
strate that the S and S" states are spin-polarized in al- 
most the opposite direction, in agreement with recent 
spin- resolved ARPES measurements [2^. The negligible 



spin-polarization around F is consistent with the overlap 
of the surface bands with the bulk projection (see Fig. [2]). 
In this area, the electron states become resonances with 
a large penetration, thus any surface effect such as the 
enhancement of the SO interaction is almost completely 
absent. 

The anisotropic character of the SO interaction is ev- 
idenced by the highly noncollinear structure of the cal- 
culated spin-polarization for S and S' . On one hand, 
we observe that the spin of each surface state is mainly 
polarized along the surface perpendicular direction, a 
phenomenon that extends to a significant area around 
the high symmetry points M and M'. Such an impor- 
tant contribution of the out-of-plane magnetization is 
a consequence of the strong in-plane gradients of the 
ionic potential, as reported in the Tl/Si(lll) surface 
[l9| . On the other hand, our calculations further identify 
an important area of almost pure in-plane circular spin- 
polarization for each state around high symmetry point 
K. 

Fig. [3j; shows the results for the calculated spin-flip 
transition probability Pss'(k) = lA/55/ (k)p/| AqP for 
an R-circularly polarized external field, Aq = Aq{x + 
iy)/\/2. Results for other polarizations, although being 
slightly different, do not substantially modify our general 
conclusions. We deduce from Fig. [Sj: that the spin-flip 
transition probability is negligible around the F point. 
This is consistent with Eg. [4] since the surface states 
are spin-degenerate here, [27| therefore SO driven effects 
such as the spin-flip excitations are weak. On the other 
hand, the most important message of Fig. [S}: is the ex- 
treme localization of the spin-flip transition probability 
in hot spots close to the high symmetry point K. A 
careful comparison of Figs. [3^ and [3)3 with Fig. [3t re- 
veals that the hot spots are localized in the boundaries 
separating the surface-perpendicular and surface-parallel 



spin-polarized regions in the Brillouin zone. This feature 
is fully consistent with Eq. |4] since the velocity operator 
introduces a k-space derivative [ii, 27| which effectively 
measures the variation of the entire wavcfunction and the 
spin-polarization through Eq. [5l Noteworthy, this sin- 
gular effect is completely absent in Rashba-like systems 
exhibiting a smooth behavior of the spin-polarization P- 

III- _ 

The light absorption rate associated to surface spin-flip 
excitations is given by the following expression, 
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w 



(6) 



Above, 7ss'(w) is the spin-flip transition rate (Eq. [3]), 
Hu! the energy of the external field and W = c\Eo\'^/8tt, 
the incident optical power per unit area. 

Fig. S] illustrates the calculated absorption rate 
associated to the spin-spit surface states of the 
A/3Pb/Ge(f If ) and the Rashba-like Au(fff) surfaces. 
For ■\/3Pb/Ge(f f f ), the absorption spectrum is bounded 
in the 0.1-0.3 eV energy range, corresponding approxi- 
mately to the spin-splitting at the Fermi level. At this 
surface, the spectrum presents a prominent peak close 
to 0.17 eV, where the spin-flip absorption rate reaches a 
remarkable maximum value of 6%. 

This result demonstrates that a significant part of the 
incoming radiation is dissipated exclusively in the spin- 
flip phenomena. The SO coupling making the orbital and 
spin degrees of freedom interrelated causes a question 
on the angular momentum conservation. As the absorp- 
tion of circularly polarized light induces a net transfer of 
angular momentum, the resulting magnetization of the 
surface states can produce a corresponding angula mo- 
mentum reservoir in a spin-orbit coupled system. Due 
to the anisotropy of the spin-polarization (see Fig. 3), 
the absorption spectrum exhibits a substantial variation 
as a function of the light polarization. As an example, 
close to 0.25 eV the spin-flip absorption rate for L po- 
larized light is approximately three times stronger than 
for R polarized light. In contrast, we find that the ab- 
sorption spectrum of Au(lll) is practically polarization 
independent. The reason is that this system exhibits an 
almost perfect isotropy, as assumed in the Rashba model. 
The magnitude of the spin-flip contribution for Au(lll) 
is relatively weak (0.8% maximum). Thus, the spin-flip 
absorption in -\/3Pb/Ge(lll) is one order of magnitude 
stronger in comparison to Au(lll). 

It is noteworthy that the bare spin-flip contribution 
to the absorption rate in •\/3Pb/Gc(lll) is around three 
times stronger than the total absorption of a graphene 
layer (2.3%), where the electron spin does not play any 



signiflcant role [341 135| . Thus, the a priori weaker rel- 



ativistic SO interaction in •y3Pb/Ge(lll) exceeds the 
effect of the usually predominant non-relativistic terms 
such as the electric dipole mechanism. The reason why 
the spin-flip contribution in this system is so important 
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FIG. 4: (color online) Calculated spin-flip absorption rate in 
\/3Pb/Ge(lll) and Au(lll). Solid (black), dashed (green), 
dotted (red) and dashed-doted (blue) lines represent the re- 
sults corresponding to the R and L circularly polarized and 
X and y linearly polarized light, respectively. Note that 
the absorption rate for -\/3Pb/Ge(lll) depends on the light 
polarization, while in Au(lll) it is practically polarization- 
independent. 



is that the hot spot matrix elements (Fig. ^) lie just in- 
side the Brillouin zone area where the S state is occupied 
and S' state remains empty. In this way, the Fermi occu- 
pation factors allow electron transitions precisely where 
the matrix elements are maximal. From the discussion 
above we can conclude that a large anomalous feature as- 
sociated to the enhanced spin-fiip excitation mechanism 
should be accessible by infrared optical spectroscopy in 
the 0.1 — 0.3 eV energy range. 

In summary, we investigate the role of the spin- 
orbit interaction on the light absorption rate of the 
A/3Pb/Ge(lll) surface. Our calculations incorporate the 
full spinor wave function structure from first-principles, 
making use of a precise integration procedure through 
a Wannier interpolation scheme for the spin-flip matrix 
elements. We flnd that a substantial part of the low- 
energy absorption spectrum is dominated exclusively by 
the spin-flip excitations associated to the spin-polarized 
surface states. Noteworthy, these transitions capture as 
much as 6% of the total incident power, representing an 
enhancement of one order of magnitude in comparison 
to the Rashba-like prototype Au(lll) surface. The ori- 
gin of such a huge absorption rate is closely related to 
the strong anisotropy exhibited by the spin-polarization 
structure connected to the surface states. As demon- 
strated in this Letter, a clear fingerprint of the spin-flip 
absorption mechanism should be accessible in the optical 
range, posing a challenge for further experimental work. 
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